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Abstract. Microquasars are potential candidates to produce a non-negligible fraction of the observed galactic 
cosmic rays. The protons accelerated at the jet termination shock interact with the interstellar medium and may 
produce detectable fluxes of extended emission at different energy bands: high-energy and very high-energy 7-rays 
produced by neutral pion-decay, synchrotron and bremsstrahlung emission in a wide energy range generated by 
the secondary electrons produced by charged pion-decay. We discuss the association between this scenario and 
some of the unidentified EGRET sources in the galactic plane. 
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1. Introduction 

The high-energy 7-ray instrument EGRET detected about 
170 sources that still remain unidentified. Among them, 
an important fraction are non-variable galactic sources. 
Several astrophysical objects could explain the nature of 
this set of steady galactic 7-ray sources: pulsars, supernova 
remnants, molecular clouds, etc. (see Torres et al. |2003. 
and also Nolan et al. 2003). However, new types of galactic 
sources could be also the counterparts of this 7-ray steady 
emission detected by EGRET. In particular, microquasars 
(MQs) have turned out to be likely sources of 7-rays 
(Paredes et al. 12000(1 . and they could be indirect sources of 
non-variable (at EGRET timescales) 7-ray emission too. 
Indeed, these objects present powerful jets that show per- 
sistent or transient activity and can carry huge amounts 
of kinetic energy during the life of the MQ (millions of 
years for high-mass MQs and longer lifetimes for low-mass 
ones, Tauris & Van den Heuvel 2003;. Following the stan- 
dard model for jet formation in MQs (i.e. Meier J2003), 
the jets of microquasars transport a hadronic component 
extracted from the accretion disk, which is formed by the 
matter accreted from the companion star. Once the jet 
terminates, two different scenarios are possible: a smooth 
transition jet-interstellar medium (ISM), where jet pro- 
tons would be released in the ISM with velocities simi- 
lar to the original jet velocity, or a shock-like interaction 
with the ISM, where jet protons would be accelerated 
by the Fermi mechanism. If the first possibility occurs, 
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these particles could be a minor but non-negligible part 
of the cosmic rays produced in our Galaxy at low ener- 
gies (about GeV energies). Moreover, in the second case, 
protons could reach higher energies, depending on the con- 
ditions within the jet-ISM shock region (for an extended 
discussion, see Heinz & Sunvaev I2QQ21. From electrody- 
namical arguments (see, i.e, Aharonian et al. l2002[) . if the 
magnetic field and the size of the final part of the jet is 
similar to that, for instance, considered in Paredes et al. 
( 2002 ) , proton energies of about 10 5 GeV may be achieved. 
Assuming that a non-thermal population of high energy 
protons is released from the terminal part of a jet, these 
particles will diffuse through the interstellar medium. At 
certain distances from the accelerator, the protons can in- 
teract with high-density regions (i.e. molecular clouds), 
producing extended and steady emission at different ener- 
gies. Due to propagation effects, the outcomes from such 
interactions can differ strongly depending on the age, 
the nature (impulsive or continuous) of the accelerator 
and the distance between the accelerator and the cloud 
(see Aharonian & Atovan H996f) . It should be noted that 
hadronic 7-rays could be produced at much smaller scales 
as well, e.g. at interaction of the relativistic jet of the com- 
pact object with the dense wind of the stellar companion, 
as suggested by Romero et al. (2003 ). 

In this work, we explore the possibility that MQs could 
initiate indirect and steady sources of 7-rays and lower 
frequency radiation throughout interactions between high 
energy protons released from their jets and nearby molec- 
ular clouds. We have developed a time-dependent model 
that calculates the broadband spectrum of the emission 
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coming out from the pp primary interactions as well as 
the emission produced by the secondary particles (elec- 
trons and positrons) created during the first process. This 
model takes into account the propagation effects due to 
both energy-dependent diffusion and the pp interaction 
energy losses of protons during their travel from the ac- 
celerator to the cloud. From this point, when referring to 
an impulsive MQ, we mean a transient jet of MQ inter- 
acting with a cloud, and by a continuous MQ, we mean 
the same but for a persistent jet of a MQ. 



2. MQ-cloud interactions and production of 7-rays 

Jets of MQs must end somewhere, although the details 
of this are still unclear (Heinz 2002). Some MQs present 
persistent jets that emit at different energy bands up to 
hundreds of AU with typical jet kinetic energy luminosi- 
ties from 10 37 to 10 39 erg s" 1 (LS 5039, Ribo 2002, SS 433, 
Marshall et al. EE 



2J Spencer Il984[) . Other members be- 
longing to this type of galactic jet sources present very 
powerful transient ejections that can extend larger dis- 
tances and contain several orders of magnitude more jet ki- 
netic energy than the persistent cases (i.e. GRS 1915+105, 
Mirabel & Rodriguez lTOOl Cygnus X-3, Marti et al. 12000")) . 
Finally, some MQs present extended structures at radio 
frequencies (IE 1740.7-2942, Mirabel & Rodriguez IT09"9")l 
or at X-rays (XTE J1550-564, Corbel et al. 120021 Kaaret 
et al. I2003fl at scales of 0.1-1 pc, which could be related 
to acceleration of particles of a jet interacting with the 
environment (i.e. Wang et al. |2003). 

Below, we initially assume that the jet has a signifi- 
cant population of protons. These protons can be cold, or 
cooled during their propagation through the jet due to, 
i.e., adiabatic losses. They can be cold in the reference 
frame of the jet, but in the frame of the ISM their en- 
ergy could be as large as Tm p c 2 , where T is the Lorentz 
factor of the jet and m p is the proton mass. However, 
they can also be accelerated in the terminal part of the 
jet owing to shocks between the matter of the jet and 
the ISM. Therefore, a population of non-thermal protons 
can extend up to very high energies, and these particles 
are released since the low magnetic field strength is not 
enough to keep them within the accelerator. Once the 
particles break free, they diffuse in the ISM with a dif- 
fusion speed depending on their energy and, close to or 
further from the releasing point, the protons interact with 
higher density regions. The interaction between the high 
energy protons and the interstellar hydrogen nuclei will 
produce charged and neutral pions (n~ , 7r + and 7r°); the 
first will decay to electrons and positrons and the second 
to photons. The primary radiation, 7r°-decay photons, is in 
the 7-ray band, but the secondary particles can produce 
significant fluxes of synchrotron (from radio frequencies 
to X-rays) and bremsstrahlung emission (from soft 7-rays 
to TeV range), and generally with less efficiency, inverse 
Compton (IC) emission through interaction with the am- 
bient infrared (IR) photons. 



2.1. Evolution in time and distance of the proton 
energy distribution 

To calculate the emissivity of a source due to different 
mechanisms, we need to know first what is the energy 
distribution of the high energy protons at different times 
and distances from the acceleration site. Therefore, we 
have assumed that the (initial) energy distribution of the 
accelerated protons follows a power-law plus a high-energy 
cutoff: 



f p0 (E p ) = KE p PeM-E P /E pmax ), 



(1) 



where K is the normalization constant of the function (to 
compare with the total number of particles or total en- 
ergy contained by them), E p is the proton energy and 
-Epmax is the energy where the exponentical cut-off starts 
to dominate over the power-law. -E pmax is determined by 
the efficiency of the accelerator and will be treated here 
as a parameter, avoiding detailed treatments of the accel- 
eration region. 

To determine the energy distribution of protons in dis- 
tance and time (f p (E p , R,t)), we have used the solution 
found by Atoyan et al. ([1995) to the diffusion equation 
(Ginzburg & Syrovatskii 1964) in the spherically symmet- 
ric case: 



dfp 
dt 



D d 
R^dR 
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dR 
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dE, 



(Pf) + Q- 



(2) 



Here, t is the source age, R the distance to the source 
of protons, D is the diffusion coefficient, which is as- 
sumed as D(E p ) oc Eg, P is the continuous energy loss 
rate (P = —dE p /dt) and Q is the source function, such 
that Q(E p ,R,t) = f p0 (E p )S(R)S(t). <5(R) and 6{t) are 
delta functions in space and time, respectively. To ac- 
count for the energy losses during the diffusive trans- 
port of protons, since the dominant cause of loss is nu- 
clear interactions, we have P{E p ) = Ep/T pp , where r pp « 
6 x 10 7 (n/l cm -3 )" 1 yr, and n is the hydrogen density 
of the medium. We assume an acceleration spectrum con- 
stant in time given by Eq. ljT]l. 

In our particular case of a power-law plus a high-energy 
cutoff injection spectrum, and a diffusion coefficient that 
depends on energy following a power-law, the general so- 
lution is reduced to: 



fp(E p ,R,t) 



KE-PeM-E P /Ep ma , x ) 



x exp 



(p-l)f 



'pp 



R 2 
R 2 



(3) 



where i?dif is the diffusion radius and corresponds to the 
radius of the sphere up to which the particles of energy 
Ep effectively propagate during the time t after they left 
the acceleration site, 

RMEp, t) = 2 (5(£ p )t^|^j . (4) 

D(Ep) = L»io(£ , p /10 GeV) x , and D w is the diffusion co- 
efficient normalization constant in the ISM. For further 
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details regarding diffusion of particles in the ISM, see 
Aharonian & Atoyan (1993). 

From Eq. we can determine the proton energy dis- 
tribution at any t and i?, starting with an impulsively in- 
jected population of particles (Eq. Q). To model the en- 
ergy distribution of protons for a continuous MQ, we will 
need to integrate in time from the initial to the present 
time, taking into account that the source function repre- 
sents the number of released particles per energy and time 
units, and not the number of released particles per energy 
unit, as in the impulsive case. 

2.2. Emission from interacting protons and hydrogen 
nuclei 

We calculate the radiative effects of diffusing high energy 
protons that penetrate a high density region of the ISM. 
We will assume that the channels of pp interaction that 
yield either 7r°, ir~ or ir + are roughly equiprobable, thus 
the functions that represent the photon, the electron and 
the positron energy distributions resulting from the pion 
decay will be taken to be similar (in the relativistic regime 
of secondary particles). 

The emissivity of the photons produced by 7r°-decay 
(g 7 ) at given R and t can be calculated through Eq. (JSJ): 



q-y (Ej , R, t) — 2 



tmax qir (E v ,R,t) 
E„ mm y/E* - mid* 



dEn, 



(5) 



where E n and E„ are the energy of the emitted photon and 
the decaying pion respectively, E^mm = E^ + m\c^ /AE 1 , 
m-x is the pion mass and c is the speed of light. The neu- 
tral pion emissivity from pp interaction (q n ) can be found 
through Eq. ©: 



q w (E w ,R,t) = cn H 



5(E„ - KE kin ) 



xa pp (E p )n p f p (E pi R, t)dE p 
o-pp(m p c 2 + £^ 7r /fc 7 r)/p(m p c 2 + E^/k^,R, t), (6) 



where k n is the mean fraction of the kinetic energy (i?kin = 
E p ~m p c 2 ) of the proton transferred to the secondary pion 
per collision, tih is the density of particles of the high den- 
sity ISM region, and a pp (E p ) is the cross section of the pp 
interaction. This approach of emissivity calculation of n - 
decay 7-rays gives quite good accuracy down to energies 
~1 GeV (see Aharonian & Atovan I2UU0"|| . 

Once we have obtained the emissivity, we need only to 
multiply by the total volume of the interaction region (V) 
to find the total number of photons generated per energy 
unit and second. We suppose that the cloud is homoge- 
neous. Also, its size must be small enough compared to R 
to assume the same energy distribution for the high energy 
protons everywhere in the interaction region. Actually, one 
can assume that for a high density region with size smaller 
than (but similar to) R, the changes in the proton spec- 
trum are not significant enough to be taken into account 
at our level of approximation. Moreover, for the involved 



timescales, protons do not interact several times within 
the cloud, because the pp interaction cooling time is longer 
than the escape time in the proton energy range relevant 
here. Then, the Spectral Energy Distribution (SED) of the 
emitted 7-rays due to the pp interaction within a cloud at 
a given t (age of the proton source) and R (distance from 
the proton source to the cloud) is: 



E 1 L(E 1 ) = E*Vq y (E y ,R,t). 



(7) 



The number of electrons per time and energy units pro- 
duced within a cloud at a certain R and t (Q e (E e , R,t)) 
presents a similar distribution, where E e is the electron 
energy. 

2.3. Secondary particles and their emission 

We have introduced bremsstrahlung, synchrotron and IC 
energy losses for the secondary leptons in order to get 
the proper electron energy distribution. Also, we have as- 
sumed that leptons remain within the cloud after their 
creation. This is a reasonable assumption as a first order 
approach for electrons emitting via synchrotron process 
at frequencies above 10 GHz, taking into account that the 
cooling time of those electrons within the cloud is shorter 
than their escape time, assuming a diffusion for electrons 
inside similar to the one for protons outside. We note also 
that for sources younger than ~ 10 4 yr, secondaries emit- 
ting roughly below 10 GHz have not had enough time to 
escape yet. We will not study emission at lower energies, 
since this requires a more complex approximation, which 
will be performed in future work. We have determined 
the dependence on time of the electron energy adopting 
quite reasonable assumptions: homogeneous and isotropic 
velocity distribution of electrons, a homogeneous distri- 
bution of hydrogen nuclei, and constant magnetic and ra- 
diation fields (B and U respectively). Regarding radia- 
tion fields, we have taken into account the 2.7 K photon 
galactic background plus an additional 50 K photon field 
considered to be produced inside the cloud (sec Table P). 
These IC losses, although unavoidable, are not significant 
unless very high density IR fields are present, which has 
not been excluded (see Sect. I3.7(l . The electron energy is 
given by: 



E e (t) 



a riHE c o 



(a n H + E e0 C) exp[a n H (£ - t )] - E e0 C' 



(8) 



where C = {a s B 2 +a c U), a = 8.1x10 



-16 „ - 



2.37xl0~ 3 



and a c = 3.97 x 10~ 2 in cgs units. Here, ao is related to 
the bremsstrahlung losses, and a s and a c are associated 
with the synchrotron and the IC losses (Thomson regime) 
respectively, to is the initial time, corresponding to an elec- 
tron energy E e Q. 

From E e (t) and the continuity equation we can predict 
how the injected secondary particles, following an initial 
Q c , will change their energy distribution in time. Knowing 
how an injected Q c evolves in time, wc integrate the dif- 
ferent evolved values of Q Q from the initial time to the 
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present moment. Doing so, we get the total number of 
electrons per energy unit at any t and R: N c (E Cl R, t). 

We can compute N e (E e ,R,i) using numerical calcula- 
tions for any value of E c , R and t. Knowing the previous 
function, the magnetic and radiation fields, and the den- 
sity of the cloud, one can easily calculate the emission from 
these secondary particles produced via synchrotron radi- 
ation, bremsstrahlung and IC scattering (i.e. Blumenthal 
& Gould Q]|7UJ|. 



3. Performance of the model for different cases 

We go through different cases that might be relevant. 
The adopted parameter values are presented in Tabled 
although some of them could change depending on the 
particular scenario. We have computed the emission pro- 
duced by four different radiative processes: 7r°-decay, syn- 
chrotron radiation, bremsstrahlung and IC scattering. For 
the impulsive case, we have considered a very active 
source, L k ~10 39 erg s _1 , but operating a short period of 
time, ~ 30 yr, which implies a E^ ~10 48 erg (transfered to 
the accelerated protons). We note that this value for the 
kinetic energy luminosity of protons is average during the 
impulsive phase, and it takes into account periods of less 
activity and recurrent super-Eddington accretion bursts 
within these 30 yr. For a continuous MQ, we have adopted 
moderate kinetic energy luminosities for the jet (trans- 
fered to the accelerated protons) of about 10 37 erg s _1 
(see Table HJ. 

In the case of impulsive MQ, as is mentioned in Sect. [3 
for the case of SS 433, we have a persistent jet present- 
ing a kinetic energy luminosity similar to or even more 
than L k ~ 10 39 erg s" 1 (Marshall et al. I2uP21 Spencer 
I1984fl . It is a super-Eddington accretion system (as can 
be seen from the huge amount of jet kinetic energy ), but 
the ejection appears to be persistent. It could also be the 
case that a microquasar suffering recurrent or periodic 
outbursts delivered, on average, the amount of proton ki- 
netic energy assumed here. For instance, GRS 1915+105, 
which can easily reach super-Eddington accretion rates, 
presents large outbursts like the one in 1994 (Mirabel & 
Rodriguez I1994|l . with kinetic energy luminosities close to 
10 41 erg s -1 (Gliozzi et al. 1999). Otherwise, for a continu- 
ous MQ, assuming a proton acceleration efficiency of 10% 
and a jet carrying 20% of the accretion power, an accretion 
rate of 10 -8 M© yr -1 is required. Therefore, the energetic 
assumptions of this work are not rare for microquasars. 

Regarding p, we have adopted a value of 2. Cosmic ray 
spectrum is assumed to be injected with power-law indices 
of around 2.3. However, from indirect spectral considera- 
tions of radiation emitted from the jets of MQs (i.e. their 
electron spectrum inferred from synchrotron emission), it 
seems more suitable in that case to take a slightly harder 
value. The different explored situations are presented in 
the next subsections. For each case, the relevant para- 
mater values are indicated in italic. 
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Fig. 1. Spectrum of protons released from a source at 
#=10 pc after t=100 yr (solid line) and 10000 yr (long- 
dashed line). The dotted horizontal line shows the flux of 
cosmic ray protons observed near the Earth, a) For an 
impulsive MQ. b) For a continuous MQ. 

3.1. Proton flux evolution 

R=10 pc; t=100, 10000 yr. 

We have computed the proton spectra (J(E p ) x E 2 - 75 ) 
for an impulsive and a continuous MQ together with the 
cosmic ray spectrum. The results are plotted in Fig. ^ We 
note a harder spectrum and more particles for a continu- 
ous MQ compared to an impulsive MQ at older ages. This 
is due to the continuous particle injection from the jet of 
the MQ in the first case. 

3.2. Proton and electron energy distributions 

R=10 pc: t=100, 1000 yr. 

We have computed the proton and the electron en- 
ergy distributions (J(E p / c ) x ^-V) within a cloud, plot- 
ting them together in the same graph to demonstrate their 
evolution with time. The impulsive and the continuous 
cases are shown in Figs. and |21 respectively. We note 
that, from the diffusion radius (Eq. J3J) and the time- 
dependent electron energy (Eq. (JHJ), we can recognize the 
peaks and turning points of both distribution functions at 
those energies that correspond roughly to the protons and 
electrons that dominate at that time. Another interesting 
feature of Figs. [21 and [3] is that the number of protons 
and electrons per volume and energy units in the continu- 
ous case moves closer to the number in the impulsive case 
as t is increased, since for the continuous source the pro- 
tons and electrons accumulate {accumulate for protons) 
in the cloud with a continuous proton injection along the 
MQ lifetime, whereas in the impulsive MQ protons do not 
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Parameter 


Symbol 


Value 


Diffusion coefficient normalization constant 




in 27 nrr^ 2 a' 1 

1U cm s 


Diffusion power-law index 


X 


u.o 


ISM medium density 


n 


u.i cm 


High density ISM region/cloud density 


"II 


iu cm 


Mass of the high density ISM region/cloud 


A 1 

M 


X 1U M0 


Magnetic field within the cloud 


B 


5 x 10" 4 G 


IR radiation energy density within the cloud 


U 


10 eV cm" 3 


Planckian grey body temperature (IR) 


T 


50 K 


Powerdaw index of the high energy protons 


P 


2 


Cut-off energy of the high energy protons 


■Spmax 


10 5 GeV 


Kinetic energy luminosity for accelerated protons in the continuos MQ 




10 37 erg s _1 


Kinetic energy for accelerated protons in the impulsive MQ 


E k 


10 48 erg 
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log (Kinetic energy [m a c 2 ]) 

Fig. 2. Proton and electron energy distributions (solid and 
long-dashed lines respectively) for an impulsive MQ at 
i?=10 pc. The kinetic energy axis units are m e c 2 . a) The 
source age is t=100 yr. b) The source age is 1000 yr. 



accumulate and secondary particles are injected with a 
changing rate. For longer periods of time, the impulsive 
MQ emission will eventually disappear, but the contin- 
uous MQ will be emitting at least during its whole life. 
The secondary particles accumulate also in the cloud for 
the impulsive MQ, but there is no injection of fresh par- 
ticles at the highest energies, so the spectrum becomes 
steeper than in the continuous MQ. Below a certain en- 
ergy, the electron energy distribution is completely flat. 
This comes from both the lower energy cutoff of the pro- 
ton energy distribution due to the dependence on energy 
of the diffusion in the ISM and the characteristics of the 
pp interaction and production of secondary particles. 
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- protons 

- electrons 
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log (Kinetic energy [m e c 2 ]) 



Fig. 3. The same as in Fig. [5] but for a continuous MQ 
(recall: a) t=100 yr, b) £=1000 yr). 



3.3. Overall spectrum at different ages 

R=10 pc: t=100, 1000, 10000 yr. 

We have computed the SED from radio to 7-ray ener- 
gies at different t. Both the impulsive and the continuous 
MQ cases are presented in Figs.0]andEl respectively. The 
luminosity peaks of primary and secondary particle emis- 
sion are shifted to lower energies as incoming primary par- 
ticles have less energy and secondary electrons lose their 
energy. Also, the secondaries accumulate within the cloud 
(assuming they do not escape), decreasing the difference 
between the fluxes from the 7r°-decay and bremsstrahlung 
channels. 

The comparison of Figs. 0] and [5] shows different spec- 
tral features. The spectrum for an impulsive MQ at low 
and high energies is steeper and presents lower energy 
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- Bremstrahlung emission 
pion-decay emission 
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Fig. 4. The SED for an impulsive MQ for R=10 pc at 
three different epochs: for 100 yr (curve 1), for 1000 yr 
(2), and for 10000 yr (3). The photon energy axis units in 
this and following plots are m e c 2 . 
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Fig. 5. The same as in Fig. SI but for a continuous MQ 
(recall: i=100 yr (1), 1000 yr (2), and 10000 yr (3)). 



luminosity peaks than the one for a continuous source. 
Also, the total emission increases first and decreases later 
for the first case, increasing asymptotically for the sec- 
ond case. This is due to the lack of continuous injection 
of new protons at high energies and the reduction of the 
number of protons at different energies with time for an 
impulsive MQ. In the continuous source, since protons are 
continuously arriving at the cloud following the same en- 
ergy distribution, there are cumulative effects on radia- 
tion. The previous issues make the impulsive case a good 
candidate to be detected at early stages of the source age 
(100-1000 yr) and at higher energies (i.e. detectable by the 
new generation of ground-based Cherenkov telescopes), 
and a good case to determine the energy of the progenitor 
particles of the emission (due to the weaker cumulative 
effects of the spectrum). Otherwise, the continuous case 
should be easier to detect at later stages of the source 
age (> 10000 yr) and at lower energies (i.e. detectable by 
EGRET), but the cumulative spectral effects would give 



a bias in the extrapolation of the typical particle energies 
from the peaks in the spectrum 1 . 

3.4. Spectra for different diffusive behaviors and 
proton origins 

R=10 pc (continuous MQ); steady regime; -E pma x ~ 
10 6 GeV; p = 2.75 (cosmic rays) and = 2 (continuous 
MQ). 

In this subsection, we have computed the overall spec- 
trum assuming that the relativistic protons belong to the 
sea of cosmic rays. We have computed also the spectra for 
two continuous MQ cases: an energy-dependent diffusive 
case and a energy-independent diffusive case. Since the 
cosmic ray sea can be considered to follow a homogeneous 
energy distribution in space and time, for this case we have 
assumed that there is no dependence on energy in the dif- 
fusive process, as well as an accelerator age long enough to 
reach the steady regime. In the second situation, it is pos- 
sible to compare the energy-dependent diffusion spectrum 
with the energy independent diffusion spectrum, taking 
the initial energy distribution of protons to be exactly the 
same (power-law index, kinetic energy luminosity of the 
jet, etc.). Again, the age has been taken long enough to 
reach the steady regime for relativistic protons within the 
cloud. In Fig. upper panel, the cosmic ray case is shown, 
the two continuous MQ cases being presented in the lower 
panel. Since galactic cosmic rays can reach higher energies 
than the adopted maximum energies for protons, we have 
now taken -E pma x ~ 10 6 GeV in the three SED to compare 
them properly. 

Fig. El demonstrates the importance of energy- 
dependent diffusion in obtaining the correct spectral 
shape. Without energy dependence in the diffusive pro- 
cess, the fluxes grow at all the wavelengths in such a way 
that there is neither spectral steepening nor emission peak 
shifting with time. Also, comparing the continuous MQ 
cases with the cosmic ray one, there is an important dif- 
ference in the spectral slope between the case with no en- 
ergy dependence in diffusion and the cosmic ray case, but 
not between the one with energy dependence in diffusion 
and the cosmic ray case. In the context of our model, it 
means that the origin of the cosmic ray sea is a group of 
sources which are far enough from their corresponding in- 
teracting clouds to suffer steepening in the proton energy 
distribution, following these protons initially a power-law 
of index ~2 (like the typical value adopted in this work). 

3.5. Overall spectrum at different distances 

R=10, 30, 100 pc; t = 1000 yr. 

In Figs. and |S1 we present the computed SED from 
radio to 7-ray energies of an impulsive and a continuous 
MQ, respectively. Due to the effects imposed by the prop- 
agation of particles in the ISM, there is a relationship 



1 Actually, for secondaries, what should be studied is the 
turning point in the spectrum due to energy losses. 
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Fig. 6. a) SED from the galactic cosmic rays interact- 
ing with a cloud. The power-law index for proton energy 
distribution has been taken to be 2.75, as is found near 
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from a continuous MQ: in one case the diffusion is energy 
dependent (solid line), in the other case it is not (long- 
dashed line). 
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Fig. 7. The SED for an impulsive MQ for £=1000 yr at 
three different R: for 100 pc (curve 1), for 30 pc (2), and 
for 10 pc (3). 



between time and distance (see Eq. (0J). This fact im- 
plies that Figs. and 03 will present similar features to 
Figs. 0] and like the decreasing difference in luminosity 
between 7r°-decay and bremsstrahlung emission and the 
shift of the emission peaks to lower energies. 

3.6. Low and high magnetic fields 

R=10; t = 100, 10000 yr; B=10~ 3 , 10~ 5 G. 
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Fig. 8. The same as in Fig. but for a continuous MQ 
(recall: i?=100 pc (1), 30 pc (2), and 10 pc (3)). 



We have calculated the SED of an impulsive MQ tak- 
ing a low and a high value for B. For a young source 
case, plotted in Fig. and a low B, owing to the small 
loss rate of secondary particles, the position in energy 
of the maximum of emission for the secondary particles 
mainly depends on the typical energy of the arriving pro- 
tons at that R and t. For a young source and a high B 
(same Fig. we recover the strong dependence of the 
secondary particles on time due to the losses. Otherwise, 
for an older source case, plotted in Fig. and a low B, 
the peak for bremsstrahlung shows at what energy the 
electrons are concentrated (roughly the same as protons) 
although the synchrotron emission peaks at higher ener- 
gies. For an older source and a high B (same Fig. I10|l . 
the computed SED again shows strong dependence of the 
secondary particles on time due to the losses, presenting 
a softer spectral slope than in the case of a small B at 
the same age. In all the previous situations it can be seen 
how bremsstrahlung and synchrotron are processes more 
or less important to each other depending on the syn- 
chrotron losses: for high B, synchrotron dominates, for 
low B, bremsstrahlung dominates. Of course, the mag- 
netic field does not affect 7r°-decay emission. 

In previous subsections the IC component has not been 
shown. This is due to the low emission fluxes produced via 
IC scattering between secondaries and IR photons within 
the cloud for the adopted density of such photons (see 
Table. QJ, much lower than the ones predicted for the 
other components (synchrotron, bremstrahlung and pion- 
decay). Nevertheless, in the next subsection, a study about 
the importance of the IC radiation is presented. 

3.7. On the importance of the IC radiation 

R=10; t = 1000; U=l, Iff eV cm- 3 . 

In Fig. ^2 we show the importance of the IC emis- 
sion comparing with the radiation from other components 
(synchrotron, bremsstrahlung and 7r°-decay processes) de- 
pending on the energy density of the infrared seed photon 
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field. We recall that the IR spectrum has been taken to 
be a Planckian grey body, with a temperature of about 
50 K. From Fig. ^2 it is seen that a high density IR pho- 
ton field (> 10 4 eV cm -3 ) is necessary within the cloud 
to have a relevant IC component in radiation. Thus, for 
clouds with standard IR photon energy densities, the dom- 
inant source of electron energy losses would be the mag- 
netic field, and for clouds with IR photon energy densities 
about four orders of magnitude higher than the standard 
value, the dominant source of electron energy losses would 
be IC scattering, with IC emission levels similar to those 
generated by the bremsstrahlung process in the same en- 
ergy range. However, IR field densities of 10 4 eV cm -3 are 
too high, such a situation being very unlikely. Of course, 
7r°-decay emission is not affected by the IR photon field 
density. It is interesting to note the small decrease of syn- 
chrotron and bremsstrahlung emission due to the moder- 
ate increase of losses because of the higher radiation field 
density, compared to the strong increase of IC emission. 
It seems very unlikely that IC losses could be important 
enough to control by themselves the radiative processes 



and the evolution of secondary particles in a reasonable 
scenario. 

3.8. Continuous MQ with a strong outburst 

R=10 pc; t = 10 5 yr (continuous MQ) and t = 200 yr 
(impulsive MQ). 

Now we combine two kind of sources, an impulsive MQ 
and a continuous MQ, to obtain the complex SED that is 
presented in Fig. 1121 To find clouds presenting emission 
with different origins is quite probable. Taking into ac- 
count the long-term variable nature of microquasars (the 
short-term is not taken into account here because it does 
not affect large-scale objects like clouds), we can find dif- 
ferent periods of activity such as persistent and impul- 
sive phases. Therefore, a spectrum like the one shown in 
Fig. H2l might be quite typical. Due to obvious differences 
in evolution for both types of emission, that originating in 
the continuous phase and that originating in the impul- 
sive phase, we find that the spectral shape presents two 
different peaks and complicated slopes, being detectable 
at radio, X-ray, EGRET and TeV enegies. 

3.9. The case of two clouds 

R=5 pc and R=50 pc, with V ~ 3.4 and 3400 pc 3 , respec- 
tively. 

In Fig. 1131 we have computed the spectrum that would 
be observed if there were two clouds at different known 
distances from the origin of the particles. We assume an 
impulsive MQ as being more clear than the continuous 
one to determine its physical parameter values due to the 
reduced particle accumulation effects. All the other pa- 
rameters except the distance from the clouds to the MQ 
would be unknown a priori, and our task is to estimate 
them from the observed spectrum using our model. 

From the characteristic frequency of the synchrotron 
emission (v c = 6.3 x lO 18 ^^) and the maximum en- 
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Fig. 12. SED of a complex source: a continuous MQ 
(long-dashed line) plus an impulsive phase (dotted line). 
i=10 5 yr for the first one and t=200 yr for the second 
one, with R=10 pc. The added spectrum is presented as 
a thick solid line. The kinetic energy luminosity for the 
continuous MQ has been taken to be 10 37 erg s , and 
the kinetic energy for the impulsive MQ has been taken 
to be 10 48 erg. 

ergy reached by both the synchrotron and the 7r°-decay 
emissions, we can roughly determine the magnetic field 
within the clouds (5 x 10~ 4 G) and the electron energy at 
the peaks and turning points. In addition, Eqs. 0]and|Bl 
can help to estimate the time of travel for the first par- 
ticles from the MQ to the cloud and the time during 
which the secondary particles have evolved (the total age, 
£=1000 yr). From previous parameters, we estimate the 
diffusion coefficient constant (Diq = 10 27 cm 2 s _1 ). The 
spectral shape can also give information about the proton 
energy distribution. Knowing the previous set of parame- 
ter values and the volume of the cloud, it could be possible 
to find the remaining parameters: riH and the total kinetic 
energy of the released protons. 

4. MQ, clouds and non variable way sources 

Statistical studies of the EGRET sources in the galac- 
tic plane point to the association between these sources 
and high density regions in the inner spiral arms includ- 
ing star formation regions and molecular clouds (Romero 
et al. 119991 Bhattacharya et al. l2003|l . It has been claimed 
that the low latitude EGRET sources can be divided in 
two different subgroups, one formed by those sources that 
show variability, and another formed by dubious or steady 
cases (Romero et al. 2004). High-mass MQs have been 
proposed as possible variable EGRET source counterparts 
(Kaufman Bernado et al. 2002, Bosch- Ramon et al. 2004). 
The results of this paper show that MQs could be indi- 
rectly responsible for some of the non-variable EGRET 
sources. Taking into account that our model can repro- 
duce a typical EGRET spectrum (see Fig. lT4*|) of a steady 
source, and the association between high-mass MQs and 
high density regions, the scenario contemplated here could 



be the explanation of, at least, a significant fraction of 
the non variable high-energy 7-ray detections. In addi- 
tion, these sources could produce detectable emission at 
TeV energies. 

In Fig. El we show the predicted SED produced by 
the interaction of protons released from a persistent jet of 
a MQ and hydrogen nuclei of a giant molecular cloud of 
10 5 Mq at R=10 pc and t=10 6 yr. At a typical Galactic 
distance, say 4 kpc, we obtain computed levels of emis- 
sion in the EGRET range similar to the ones observed by 
this instrument, with spectral slopes similar to the average 
value (~2.16) of this steady subgroup of EGRET sources 
of the galactic plane (see Bosch- Ramon et al. l2004fl . These 
points make it reasonable to suggest an indirect associa- 
tion between the steady EGRET population and MQs. We 
note that, due to the small distance between the cloud and 
the MQ considered here in comparison to the cloud size, 
our results are a first order approach and an inhomoge- 
neous model is required to improve the accuracy of our 
predictions, which will be presented in a future work. 

MQs can emit 7-rays by themselves (i.e. Bosch- Ramon 
et al. I2004|l . The 7-ray emitting MQ and cloud could have 
an angular separation too small to be distinguished by 
EGRET (although they could be distinguishable at lower 
energies). In such a case, the stronger emitter of 7-rays 
would dominate the spectrum and variability properties 
and the next generation 7-ray instruments, with better 
angular resolution than EGRET, would be required to 
separate them. 

One interesting issue is the shape of the young proton 
spectra (see Fig. ^) and the subsequent electron energy 
distributions (see Figs. El and I3J). As is seen from the men- 
tioned plots, there is a natural low-energy cutoff in the 
proton spectrum produced by the energy-dependent diffu- 
sion, which implies an almost monoenergetic distribution 
of electrons peaking at high energies. This characteristic of 
the secondary particle energy distribution has a clear im- 
plication for the related synchrotron spectrum: a slope like 
the one presented by a monoenergetic electron spectrum 
at low frequencies with a spectral index of 1/3. Therefore, 
in the context of a MQ interacting with a cloud, we expect 
to observe extended and steady emission with a spectral 
index of 1 /3 at radio and even shorter wavelengths. For in- 
stance, detected radiation coming from the region near the 
galactic center shows a similar spectrum (see, i.e., Lesch 
et al. I1988J) . possibly generated in a scenario like the one 
proposed here. 

5. Conclusions 

The study of the emission coming from high density re- 
gions of the ISM provides us with information about accel- 
erators of high energy particles located near these clouds. 
In this work, we have explored the implications of the 
presence of a MQ close to a cloud. The obtained results 
predict that extended radiation from radio frequencies to 
soft X-rays, generated by synchrotron radiation, could be 
detected from clouds, together with GeV and TeV 7-rays 
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produced mainly via 7r°-decay. Further information can 
be obtained about cloud internal characteristics (average 
magnetic fields, average densities). If the jets of MQs are 
sources of relativistic hadrons, then nearby clouds can be 
indicators of such hadronic accelerators with a specific 
broadband spectrum from radio to very-high energy 7- 
rays. 
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Fig. 13. SED of the observed spectrum from two clouds 
interacting with an impulsive MQ at distances 5 and 50 pc 
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Fig. 14. SED of a continuous MQ at t=l0 6 yr and 
R=10 pc reproducing a typical GeV source (solid line). 
The value of the cloud mass is 10 5 M©. An extrapola- 
tion of a typical EGRET spectrum up to ~10 GeV is also 
shown (long-dashed line). 



